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Abstract. In this study, we propose a recursive approach to study the transport properties of atomic 
wires. It is based upon a real-space block-recursion technique with Landauer's formula being used to 
express the conductance as a scattering problem. To illustrate the method, we have applied it on a model 
system described by a single band tight-binding Hamiltonian. Results of our calculation therefore may 
be compared with the reported results on Na-atom wire. Upon tuning the tight-binding parameters, we 
can distinctly identify the controlling parameters responsible to decide the width as well as the phase of 
odd-even oscillations in the conductance. 
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1 Introduction 

A mono-atomic quantum wire has a cross-section of one 
atom and is several atoms long. Such a system can be 
formed by pulling atomic contacts using a scanning tun- 
neling microscope (STM) or a mechanically controllable 
break junction (MCBJ). The experimental evidences of 
formation of such atomic wires have been reported by 
Yanson et al .L and Ohnishi et al [2]. However, issues 
involving electronic structure, transport properties and 
the influence of contacts with macroscopic conductors are 
still not fully settled. The most striking feature of such 
a mono-atomic wire is the non-monotonic behavior of its 
conductance as a function of the number of atoms along its 
length. Oscillatory behavior of conductance as a function 
of the length of the wire has been experimentally observed 
in Al, Pt and Ir wires by Smit et al [3] , in Au and Ag wires 
by Thijssen et al [I], in Na wires by Krans et al [5] and 
in wires of atoms of higher valency by Yanson et al pQ. 
This oscillatory behavior has been attributed to interfer- 
ence effect resulting from the changes in the connection 
between the wire and the edge of the electrode when new 
atoms are pulled into the wire. 

Theoretically, the conductance of atomic wires can be 
calculated using the Landauer formalism in which one 
relates the conductance G in the linear response regime 
to the transmittance at the Fermi energy, T(Ep), as : 
G = (e 2 /ttTi)T(Ef). There are several numerical meth- 
ods to calculate transmittance. They may be broadly clas- 
sified into wavefunction and Green function methods. In 
wavefunction technique, one solves for scattering wave- 
function of the system using methods like transfer ma- 



trix method [51171151131110] ■ finite difference method [TT] or 
by solving Lippman-Schwinger equation [T2"lll3j . All these 
methods use various level of approximations to describe 
the electronic structure of the system, starting from semi- 
empirical models like extended Hiickel [6j[7] to fully atom- 
istic descriptions based on the density functional theory 
(DFT) of Kohn and Sham. In most of the wave function 
methods, however, the electronic structure of the scatter- 
ing region is resolved in detail, while the leads are modeled 
by a free electron gas [8,9,12,13 . Alternatively, conduc- 
tance can be calculated by Green's function method which 
does not require the explicit calculation of the scattering 
wave function. The main effort in this approach is to calcu- 
late the Green function of the central wire in the presence 
of the coupling to the leads. The effect of coupling to the 
leads is taken into account through self-energy terms. Var- 
ious methods based on this approach, mainly differ in the 
choice of basis set used to represent the Hamiltonian and 
self-energy matrix, e.g Gaussians [14j . numerical atomic 
orbitals [T3], wavelets [TB] and plane wave basis [17] . Al- 
though implementations within the Green function as well 
as the wavefunction approaches have been carried forward 
using various different techniques, it can be shown that the 
approaches are completely equivalent for non-interacting 
electrons [18]. In this paper we propose a combination of 
the scalar and vector recursion techniques [T91I201I2I] as a 
viable and efficient means to study the transport proper- 
ties within the general scheme of wavefunction approach 
of a lead-wire-lead system. 

As a first step, starting from the junction point of the 
leads and the wire, we have used the recursion method of 
Haydock et al [19, 2D] to map the quasi-lD semi-infinite 
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leads with finite extensions along lateral directions onto 
equivalent chains. In the process, the information of the 
shape of the leads is encoded in the recursion coefficients. 
The recursion coefficients converge and the asymptotic 
part of the chains (or terminators) resembles periodic one- 
dimensional leads. This naturally divides the whole system 
unambiguously into an effective scattering region and two 
attached ballistic, periodic chain leads. The initial varia- 
tion of the recursion coefficients contributes to the scat- 
tering region. 

In the second step, the scattering matrix has been cal- 
culated by solving the Schrodinger equation using the vec- 
tor or block recursion method of Godin and Haydock [21] 
and applying the wave-function matching conditions at 
the lead- wire interfaces. 

In the final step, we have applied Landauer's formula 
to express the conductance of the quantum wire. 

The advantage of the proposed method is that it is a 
real space based method. Moreover, to solve the scattering 
problem, one does not need to calculate the wave func- 
tion explicitly. By putting the boundary conditions at the 
lead- wire junctions and at the end of vector chain, one can 
directly calculate the scattering matrix. The basic inputs 
in the above described procedure are the tight-binding 
(TB) Hamiltonians for the wire and the leads. In order 
to illustrate the method we have applied it to a model 
system described by a single-band TB Hamiltonian. We 
demonstrate the validity of this approach by comparing 
the results obtained out of a simple model system with 
the reported results on monoatomic Na-wire calculated 
from first-principles. Moreover, by tuning the TB param- 
eters of the wire and the leads we have provided a detailed 
understanding of the various features of the problem be- 
yond what has been reported earlier 31,32,33,34,35,36, 
I5T1I551I55] . The proposed technique being based on recur- 
sion, relies on the sparseness of the starting Hamiltonian. 
The coupling of this method to fully atomistic DFT de- 
scription of the electronic structure, however, can be easily 
achieved in terms of Wannier function based Hamiltonian 
constructed out of DFT calculations. 



2 Method 

2.1 Converting quasi-lD lead to effective ID chain 
and detection of effective scattering region 

Each lead with finite lateral dimension of N x x N y is de- 
scribed by a Hamiltonian in a tight-binding basis {\ifJ-)} 
where i labels a site and /i a particular channel, e.g. in a 
linear combination of atomic orbital (LCAO) type formu- 
lations, this would be the angular momentum labels (£, m) 



where x are the AT* nearest neighbors of site i on the lat- 
tice. To obtain the equivalent 'chain' or a tri-diagonal rep- 
resentation of the Hamiltonian, one needs to change the 
tight-binding basis {!«/**}} to a new one : {|n^>}, obtained 
recursively. 

Taking |1^>= |l/i) as the starting state, where 1 labels 
the site at the middle of the cross-sectional edge of the 
lead where it is connected to the wire (cf. Fig. la) and \i 
is any one of the 'orbital' indices, we generate : 

\n + 1>= flw|n> -<|n» -I3£ 2 \n - 1» (2) 

and 

M _ <n|iJi cad |n> _ <n - l|Hi cad |n> 

«n|n» [«n|n»«n- l|n- 1>] 7 " 

(3) 

The equivalent Hamiltonian in this new basis is tri-diagonal 
(chain-like) : 

oo 

#iead = «»l n >< n l + ^(|f!»«n + 1| + \n + l»«n|) 

n=l 

(4) 

where the index n labels the 'atoms' of the equivalent 
linear chain for the ji channel and a% and j3!^ signify its 
on-site and the hopping terms respectively. The sequences 
{o^,/3^} converge, so that for a given error tolerance e, 
there exists an integer c such that for n > c, \a^ — a^\<e 
and \&n — Pc \ <- e - The 'terminator' approximation puts 
a n->Pn — a ciPc f° r au n> c. The initial c 'sites' of both 
the equivalent input and output chain leads, therefore, 
contribute to the effective scattering region in addition to 
the wire. The wire of length 2M atoms is described by a 
Hamiltonian in the tight-binding basis as: 

Hjf = eJJL + ' OWi + hi-i) (5) 

The opposite ends of the wire are coupled to the semi- 
infinite effective ID leads via hopping matrix element t c . 
Note that as the starting state is same in the recursion pro- 
cess, the coupling coefficient t c remains same during the 
change of tight-binding basis set. The procedure therefore 
converts the whole system into an infinite linear chain in 
which effective scattering region is of extension 2M + 2c 
sites and the rest represents the ballistic parts of the leads 
which do not participate in the scattering process. 



^-^ ^-^ - / 2.2 Calculation of the scattering matrix 

-fflcad = 2^ 2^ H ij V> V)(hH\ 

d,il Let us consider 2M+2c = IN and for convenience, rename 

K' = &Jw> ^ + C' hi+x (1) the sites as follows (cf. FiglTb) : 
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V T W,n+1 + S n >,n-1 



n 7 fi n' ,fi 



TTfJ.fi _ ~flSj Sj 

11 nn' ~ c nW™ 



with 



The subsequent members of the basis are generated 
from 

B\\$ 2 } = {H-A^x}, 
(6) B\ +l \$ n+1 } = {H - An)\$ n } - B n \$ n _i] for n>2 

(10) 

The matrix inner product is defined as the 2L x 2L 
matrix 



n < rii left ballistic lead 

ni < n < ri2 left scattering lead 
n,2 < n < n 3 scattering wire 



a c+l-n 

r". 

wire 

a n-{c+2M) n 3 ^ n — 77-4 rl ght scattering lead 
p 



a 



n > ri4 



right ballistic lead 



and 



(7) 



n < n\ 
n\ < n < n% 
n = n 2 
n 2 < n < n 3 
n = n 3 
Pn-(c+2M) n 3 <n< r H 
^ j3g n> n 4 



wire 
if* 



ri2 



c, n 3 = c + 



left ballistic lead 
left scattering lead 
left junction 
scattering wire 
right junction 
right scattering lead 
right ballistic lead 

(8) 

2M and n 4 = 2c + 2M = 



where ni = 0, 
27V. 

We then compute the scattering matrix of the above 
mentioned system using the block or vector recursion tech- 
nique introduced by Godin and Haydock [21] . The essence 
of the vector recursion technique is the block tridiagonal- 
ization of the system Hamiltonian by changing to a new 
orthogonal set of vector basis, with the restriction that 
the ballistic part of effective ID lead Hamiltonian remains 
unchanged. In this last aspect, it differs from the stan- 
dard Lanczos method [2D] . The numerical stability of this 
method [32] has been established in studying problems re- 
lated to Anderson localization and quantum percolation 
model [33] and layering transition in 2D nano-strip 24 
previously. Below we describe the method briefly. 

Let us consider, for the sake of demonstration, that we 
have two ID ballistic leads, one incoming and and another 
outgoing connected to the opposite ends of the scattering 
region at positions |1) and |2iV) (cf. Fig [lb). The leads 
and the scatterer have fi = 1,2, ... L scattering channels. 
The member of the new basis are generated by clubbing 
the input and output leads together as the vector lead and 
so is the scattering region (cf. Fig[Tfc) . The lead states are 
chosen to be 



\n, fJ> = 1) \n, n = 2) . . . n, fi = L) 
|m, fi — 1) m, fi = 1) . . . m, /i = L) 



(9) 



with in = 2N + 1 - n and n = 0, -1, -2, . . . , oo. The 
starting state within the scattering region is chosen to be 



l*i} = 



|1,M = 1> |1,M = 2}...|1,M = £> 
\2N,n=l) \2N,fi = 2) ...\2N,fi = L) 



w®i*x}=(c£) 



where the L x L matrices A, B, C and D are 



A = 

C = {2Nfi\lfi') 



B = (lfi\2Nii'} 
D = (2Nfi\2Nfi') 



If this matrix is / then the states are called orthogonal. 

It can be shown that the 2L x 2L matrices A n and B n 
are block-tridiagonal members of the matrix representa- 
tion of the Hamiltonian in the new basis: 

An = {$„] H\$ n } B n = {$ n \ H\$ n ^ x } (11) 

so that the transformed Hamiltonian matrix can be 
divided in 2L x 2L blocks, with only non-zero diagonal 
and sub-diagonal blocks. 

The wave function may be represented in this new 
basis by a set {ip n } so that \&} =J2 n tp n \ ( l , n}- These wave 
function amplitudes ip n also satisfy an equation identical 
with (fTU|) . 

In the ballistic part of the effective ID lead chain, the 
onsite and hopping terms do not vary. Therefore, electron 
potential must be periodic in this region and the solution 
to the Schrodinger's equation in the fi channel of the ballis- 
tic leads are traveling Bloch waves : A J2m expfiim^] \m). 
As the wave travels in the leads, the phase of its wave func- 
tion changes by 9^ — cos -1 [(E ~ a£)/(2/3g)], where E is 
the energy of the incoming electron. In order to have only 
propagating solutions, we fix E as real and \E — a c \ < 2f3 c . 
This sets the energy window. This is reasonable because 
eventually only the propagating modes enter in the ex- 
pression of transmission matrix elements, as evanescent 
modes do not contribute to the transmission directly. In 
the fi channel of the incoming lead, the incident and the 
reflected waves can be expressed as a sum 



exp(- 



\rn) 



The second term is the reflected wave in the fi channel 
from incident waves in fi' channels. 

In the fi channel of the output lead there is a trans- 
mitted wave from incident waves in fi' channels |25] 

A^^V" '(E) exp(-im6' i ')\m) 

m fi' 
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r w (£j) anc f ^mm are the complex reflection and 
transmission coefficients. The boundary conditions are ther 
imposed from the known solution in the leads at the junc- 
tion labeled by and 1: 



Semi-infinite 3D input lead wire of 2M atoms Semi-infinite 3D output lead 



EA S . 



(E)} 



' £ { exp[i^] V + r^' (E) expHfl"'] } 
£ , t^'(.E)exp[-i0'*'*'] 



(12) 



(13) 



The amplitude at the n-th basis may be written as 



tp n = X n tp + Y n i>\ 



(14) 



where X n and Y n satisfy the same recurrence relation 
as (fTU|) with EI replacing H and also satisfy the boundary 
conditions Xq = I and X\ = 0, while Yq — and Fi = 
/. Note that X and Y are 2L x 2L matrices and -0 n 's are 
column matrices of dimension 2L. 

This new basis terminates after v — N steps, as the 
rank of the space spanned by the original tight-binding 
basis remains unchanged after the transformation. Hence 
the recursion also terminates after v steps. This gives an 
additional boundary condition, 



X v+ iipo + Y v+1 tp-i = 2 l> 



2L- 



(15) 



If we now interchange the incoming and outgoing leads, 

we get a similar pair of equations for r'^ and i"*** , the 
transmission and reflection coefficients for a wave incident 
from the second lead. Time-reversal symmetry demands 
that t must be the same for waves of the same energy in- 
cident from either lead so that i MM = i' w . In addition, 



and 



differ only by a phase factor. Solving these 



equations for the S-matrix in the effective scattering re- 
gion, one has, 



S(E) = -(X N+1 + Y N+1 E*{9))- 1 (X N+1 + Y N+1 E{8)) 

_ ( r w ' (E) (E) \ 
~ \t^'(E) r'{E)^' I 



(16) 



where 



E{6) 



exp(i^)5 w / 






exp(i6" i )5 A1 



The various steps of our approach are shown schematically 
in Fig. [U 

The conductance is then given by Landauer formalism. 
It expresses the electronic conductance in one-dimensional 
conductor as a quantum mechanical scattering problem 
and relates to the total transmission probability of the 
electron at the Fermi level [26]. T(Ep), as 



(17) 



(a) 



Recursion 



Input lead Output lead 

. contribution to wire of contribution to 

Ballistic lead ... . . R a kt r lead 

scattering region 2M atoms scattering region »*ti»sin. itdu 



c-1 c c+1 c+2 c+2M 



C+2M+3 2N 2N+1 



(b) 



Effective scattering region - 

1 Formation of 
vector basis 




I2N>/ I2N-1> 



(c) 

Fig. 1. Recursive reduction of (a) a system of two quasi ID 
semi-infinite leads plus ID wire into (b) a system of infinite 
linear chain. The numbering of sites of the equivalent ID infi- 
nite chain is shown in (b). The directions of incident wave e m9 
to the scattering region, the reflected wave re~ mB and trans- 
mitted wave te m6 from the scattering region are shown with 
arrows, (c) Formation of vector basis by folding the infinite 
chain and clubbing the two sites together. 



We have assumed in our derivation the leads to have fi- 
nite lateral dimensions, though in actual experiments, the 
leads are three dimensional. In order to study the influ- 
ence of the extend of the lateral dimension of the leads, we 
have checked our results with increasing size of the lateral 
cross-section of the leads. Though the quantitative values 
do change with the change of the lateral dimension, the 
qualitative nature holds good in every case. In addition, 
we find for sufficiently large choice of lateral cross-section, 
the recursion coefficients (a, /?) for the leads converge to 
that of a bulk system (see for details section 3.1), indicat- 
ing for such large lateral cross-section, the leads behave as 
truly three-dimensional leads. The other approach to the 
problem could have been inclusion of periodic boundary 
condition as have been adopted in various works P^lfTolfTTl 
H51l2"Tll2"51l2l)] . If we impose periodic boundary conditions 
on the lead surfaces, only those modes which are consistent 
with the boundary conditions can travel through them. In 
an earlier paper we have shown |30j . that it is possible 
to change over from a site to a mode basis and reformu- 
late the vector recursion in the new basis. The lead is 
broken up into slices perpendicular to its length and the 
composite site label i is partitioned into two : one is the 
slice label s and the other is the position on this slice k. 
The basis {|sfc}} is then converted into a slice- mode basis 
{|si/)} and the Hamiltonian is expressed in this new ba- 
sis. Exactly as the Hamiltonian (6) mixes orbital labeled 
channels because of off-diagonal terms , in the mode 
based formalism off-diagonal terms in the corresponding 
Hamiltonian causes the outgoing wave to be of a mixed 
mode type even if the incoming wave is in a single mode. 
Otherwise, the formal methodology is identical in the two 
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Fig. 2. (Color online) Variation of recursion coefficient /3 with 
recursion steps for different cross-sectional sizes of the lead, 
with ti e ad = 2.0 and ei ea d = 0.0. The lower left inset shows 
the variation of converged values of (/3 C ) with cross-sectional 
size. 



cases. The interested reader is referred to the paper refer- 
enced above for the details. 



3 Application 



decreases which facilitates to easier computational exe- 
cution). Variation of converged value f3 (f3 c ) with lateral 
cross-sectional size of the lead is shown within the inset. 
Notice after a certain cross-sectional size (30 x 30 with 
the present choice of parameters), (3 converges to that of 
the bulk cubic system, which is 6 for choice of ei oa d = 
and tiead = 2.0. This implies that there is no further effect 
of finite cross-section on (3 C and therefore on conductance. 
In changing cross-section from 5x5 to 30x30, f3 c changes 
by 9 % (The change is mostly coming in going from 5x5 
to 11x11 and after that change is very minute). The cor- 
responding change in conductance is much smaller : 5 % 
for even numbered wire, less than 1 % for odd numbered 
wire (cf. Fig. [3j) . We have checked that qualitatively re- 
sults are same for any cross-sectional size, but if one is 
concerned about the quantitative value, then convergence 
of /3 C with cross-sectional size has to be checked. For our 
chosen case, f3 n converged after around 80 recursions. So 
in our calculation, c was taken to be 80 and /?„ = /3so for 
n > 80. For a model simple cubic lead, the Fermi energy 
is at E l p ad = ei ca d which is obvious from the dispersion 
relation E — ei ca d + 6£i ca d cos(fca) for a half-filled band 
with Fermi wave number kp = 7r/2a. Fermi energy Ep of 
the lead-wire-lead system tends to align with the Fermi 
energy of the semi-infinite lead. Conductance is therefore 
calculated at Ep = ei ca d- 



3.1 The model 

Since our aim here is to propose a method rather than to 
evaluate the properties of any particular system in quan- 
titative detail, we present the study of a simple model 
system. It consists of a wire, few atoms long, sandwiched 
between two identical semi-infinite leads. Both the wire 
and the leads have only a single channel (L — 1) corre- 
sponding to single s-band. The cross-sectional size of each 
lead in our calculation was 5x5. The atoms in the leads 
form a simple cubic lattice. Each lead is described by a 
single-band nearest-neighbor TB Hamiltonian. Two TB 
parameters in equation (|TJ) : the on-site and the hopping 
terms, which are material specific in realistic cases, have 
been chosen as e l = ei ca d = and t^ = ii oa d = 2 in some ar- 
bitrary energy unit. As a first step, we performed scalar re- 
cursion on the leads to convert them into equivalent chains 
and determined the contribution of each lead to the effec- 
tive scattering region. We took the starting state of recur- 
sion to be the sites which bound the leads to the wire. For 
a simple cubic lead, the recursion coefficients a n = ei ea d = 
while the other coefficients j3 n fluctuate with recursion 
steps converging to an asymptotic value /3 C . The conver- 
gence of j3 n depends on the lattice structure and lateral 
cross-sectional size of lead. We have not used any peri- 
odic boundary conditions in the lateral direction. To see 
the finite size effect of lead in the lateral direction, we 
have repeated our calculations for different cross-sectional 
sizes. Convergence of /3 with recursion steps for different 
cross-sectional sizes is shown in Fig|3] As the cross-section 
increases, the number of recursion steps required to con- 
verge (3 decreases (i.e size of effective scattering region 



3.2 Results and discussions 

Our calculated conductance as a function of number of 
atoms in the wire is shown in FigSl The inter-atomic hop- 
ping within the wire was fixed to t w - n - c = 2 = fi ca d • To 
consider the effect of charge transfer between the leads and 
the wire [Ae = ei oa d - Cwhc ^ 0], we considered the cases 
\Ae\ = 0.2 and \Ae\ = 0.4 along with no charge transfer 
condition \Ae\ =0. Non zero values of Ae were achieved by 
changing e w i ro . For each case, we considered three types of 
coupling between lead and atomic wire, given by lead-wire 
hopping coefficient t c : (i) t c = 0.9 U ea d for weak coupling, 
(ii) t c = ti ead for direct coupling and (iii) t c = 1.1 U ead for 
strong coupling. 

The odd-even oscillation in conductance for no charge 
transfer situation is obvious from FigHk- i.e wires with 
odd-number of atoms have larger conductance as com- 
pared to the even numbered ones, for all the three types 
of couplings. In this ideal case of no charge transfer situ- 
ation, all odd numbered atomic wires have conductances 
equal to quantum unit Go = {e 2 /i:K) for all the three 
coupling cases, while for even numbered wires the con- 
ductances are lower. The conductance of even-numbered 
wires decreases as one moves from strong to direct to weak 
coupling. This means that the amplitude of conductance 
oscillation is a maximum for the weak coupling case (t c 
= 0.9 ti ea d) and it is a minimum for the strong coupling 
case (t c = 1.1 tiead)- Similar results have been predicted 
by Khomyakov et al for a system of ID wire connected 
with ID leads using tight-binding calculation [35]. The 
dependence of conductance on t c for even numbered wire 
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Fig. 3. Variation of conductance with cross-sectional sizes for 
an even numbered wire, N — 4 considering direct coupling 
and \Ae\ = 0.2. Variation in conductance is around 5% for 
even numbered wire. For odd numbered wire (not shown in 
the above plot), it is less than 1%. 
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It follows the same trend as observed in the top panel 
of FigQJ Another point to notice from FigHJi is that the 
length of the wire does not have any effect in this case 
either on amplitude or phase of the oscillation. All odd 
numbered atomic wires and also all even numbered atomic 
wires have equal conductances for a fixed coupling type. 
The odd-even oscillations for our single band model sys- 
tem is in agreement with the results on Na-atom wires [3T1 
!521I^IMll55pgp7j. 

For \Ae\ = 0.2 in FigQJ), the odd-even oscillation in 
conductance is again obvious within the range of our plot. 
However, a close view to these curves, indicates that the 
conductance of the odd numbered atomic wires gradually 
decreases while that of the even numbered atomic wires 
slowly increases. If we increase the number of atoms in the 
wire to more than 10 (not shown in FigH]), we find that the 
parity of conductance oscillation changes from odd-even 
to even-odd, i.e even numbered atomic wires now have 
larger conductance than the conductance of odd numbered 
atomic wires. 

Increasing |Z\e| to 0.4, there will be larger charge trans- 
fer between the lead and wire. This causes a change in the 
parity of conductance oscillation from odd-even to even- 
odd for even shorter lengths : around wire lengths of 7 
atoms as shown in Fig|4]:. From our study, it is there- 
fore clear that the amplitude of conductance oscillation 
is mainly controlled by coupling coefficient between the 
leads and the wire, while the parity of the oscillation is 
controlled by the length of the wire and onsite energy dif- 
ference between wire and lead Hamiltonians. There are 
few reports [361138] which indicate the change of parity of 
conductance oscillation and variation of oscillation ampli- 
tude for Na-atom wire. 

To understand the odd-even oscillation in the conduc- 
tance, we have calculated the density of states (DOS) and 
transmittance as a function of energy for wires of lengths 
N = 1, N = 2 and N = 3 for the case of zero charge 
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Fig. 4. (Color online) Conductance as a function of the number 
of atoms in the wire. For given values of £i ca d and t w h-o, we 
consider three cases: (a) \Ae\ = 0, (b) |Zie| = 0.2 and (c) \Ae\ = 
0.4. For each case, we consider three types of coupling between 
lead and wire : green curve for strong coupling, black for direct 
coupling and red curve for weak coupling. 



transfer. The results shown in FigJS]bring out the essential 
mechanism which has been discussed also in ref. [35] . For 
odd-numbered atomic wires, the Fermi-energy falls in the 
non-bonding peak and it gives a maximum in the trans- 
mittance at E — Ep. On the other hand, for even num- 
bered wire, the Fermi level lies in the minimum between 
the bonding and non-bonding peaks and consequently it 
exhibits a less transmission at E = Ep. 

In order to explain the flipping of conductance oscilla- 
tion from odd-even to even-odd, we have plotted in Fig|6] 
the DOS around Fermi energy of several even numbered 
atomic wires for two cases of no charge transfer and fi- 
nite charge transfer situations. When there is no charge 
transfer (|Z\e| = 0), all the even numbered wires have a 
minimum in DOS at Ep. When we allow sufficient charge 
transfer to take place between wire and lead, DOS at Ep 
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Fig. 5. (Color online) Density of states (left panel) and total 
transmittance (right panel) in no charge transfer (\Ae\ = 0) 
and direct coupling (t c = tuad) situation for wires of lengths 
N = 1, N = 2 and N = 3 (from top to bottom). The dashed 
lines label the Fermi energy. 
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Fig. 7. (Color online) Plot of conductance as a function of 
wire size for two values of \Ae\ : circles connected by solid line 
for \Ae\ = 0.6 and squares connected by dashed line for \Ae\ 
= 0.8. Conductances for all even numbered wires and all odd 
numbered wires are shown separately by red line and green line 
respectively within the insets. Inset (a) for \ At\ = 0.6 and inset 
(b) for \Ae\ = 0.8. Inset (c) shows the variation of period of 
parity flip with \Ae\ and solid line through the right triangles 
is the fitted curve. 
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Fig. 6. (Color online) DOS around Fermi energy for several 
even numbered wires in two cases - \Ae\ = (left) and \Ae\ 
= 0.4 (right). For \Ae\ = 0.4, DOS at E F gradually increases 
with wire size. For N = 14, the deep at Ef for \Ae\ = case 
is replaced by a peak at Ef for \Ae\ = 0.4. Lead-wire coupling 
was considered to be of direct type in all cases. 



gradually increases with increasing length and at a criti- 
cal size, a peak will appear in DOS at Ep. The opposite 
occurs for odd-numbered atomic wire (not shown in the 
Fig[B]) and nature of conductance oscillations flips from 
odd-even to even-odd. If we increase the wire size fur- 
ther, again after another critical size, the odd-even nature 
of oscillations is restored. This repetition of even-odd or 
odd-even oscillation continues with increasing wire size. 

To investigate the effect of |Z\e| on the period of par- 
ity flip, wc have plotted the conductance as a function of 



number of atoms in the wires in Fig(7] for two different 
values of \Ae\ keeping coupling parameter fixed. One can 
see clearly that as \Ae\ increases, frequency of parity flip 
increases i.e period of parity flip decreases. To get a better 
insight into parity flip, we have shown in Figj7]the conduc- 
tance variation of odd numbered and even numbered wires 
separately in insets (a) and (b) for \Ae\ — 0.6 and \Ae\ = 
0.8 respectively. Each curve shows sinusoidal-like varia- 
tions. Curves of even numbered and odd numbered wires 
together constitute loops. Length of one loop indicates the 
wire size required to flip the oscillation from odd-even to 
even-odd. The nodal points of loop indicate the bound- 
aries between odd-even and even-odd. Two consecutive 
loops constitute a period. Number of loops in inset (b) is 
larger than in inset (a). Variation of period of parity- flip 
with \Ae\ is shown by right triangles in inset (c). For |/ie| 
= 0, there is no parity flip i.e period is infinity. As |Z\e| in- 
creases, the period decreases. For sufficiently large value of 
\Ae\ (>> t c ), odd-even nature of conductance oscillation 
no longer persists. Period of conductance oscillation then 
changes to more than two atoms which is the characteris- 
tic of wires consisting of atoms of higher valency. Within 
the odd-even nature of conductance oscillations, the pe- 
riod of parity flip goes roughly as (A/|Z\e|) + B with A = 
12.605 and B = -0.253. 

However, coupling parameter t c does not have any ef- 
fect on period of parity flip. To check this, we have plotted 
the conductance variation with wire size in FigJ5] for two 
different coupling constant t c keeping |Z\e| fixed, while in- 
sets (a) and (b) show the conductance variation of odd 
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Fig. 8. (Color online) Plot of conductance as a function of 
wire size for two values of t c - circles connected by solid line 
for i c = 1.1 ti ea d and squares connected by dashed line for t c 
= 0.9 ti ea d- Conductances for all even numbered wires and all 
odd numbered wires are shown separately by red line and green 
line respectively. Left inset corresponds to t c — 1-1 tiead and 
right one to t c = 0.9 ti ea d ■ 



Table 1. Values of conductances for various wire lengths in 
absence of mirror symmetry. t c i and t C 2 are the two lead-wire 
couplings of the two junctions. 



wire length 


tcl — tiead 


tcl — tiead 


tcl — tiead 


(N) 


t c 2 = tiead 


tc2 = 1-1 tiead 


t c 2 = 1-2 tl ea d 


1 


0.9996 


0.9906 


0.9672 


2 


0.7164 


0.7876 


0.8482 


3 


0.9996 


0.9906 


0.9672 


4 


0.7164 


0.7876 


0.8482 


5 


0.9996 


0.9906 


0.9672 



numbered wires and even numbered wires separately. The 
number of loops in both the cases are same indicating t c 
has no effect on the period. Larger width of loops in inset 
(b) compared to that in inset (a) indicates that t c controls 
the amplitude of odd-even oscillations . 

So far we have studied the role of the charge neutrality 
on the conductance oscillations of monoatomic wires con- 
sidering the mirror symmetry between the two junctions. 
We found odd-even oscillation in the conductance. More- 
over, for no charge transfer situation, conductance of odd 
numbered wires is quantized to Go, while for even num- 
bered wires conductance is less than Go . Now we consider 
the situation where mirror symmetry between two lead- 
wire junctions is broken by using two different coupling 
coefficients for two junctions, keeping charge-neutrality in- 
tact. Table Q] contains our result. Clearly, odd-even oscil- 
lations still occur, but the conductance quantization for 
odd numbered wires is weakened by reducing the value 
less than Go- This observation is in accordance with the 
previous study (3"3"ll4*0] . 



To conclude, we have used combination of real-space 
based scalar and vector recursion techniques to study the 
transport properties of a lead- wire-lead system. Our study 
on model system described by single band TB Hamilto- 
nian provides a detail understanding of the effect of lead- 
wire coupling on the conductance of monoatomic wire. 
Working with model system, gives us the freedom of chang- 
ing the model parameters and allows us to work with much 
longer wires than usually considered in liturature. Odd- 
even oscillation in the conductance with increasing length 
of wire has been observed in agreement with earlier studies 
31 , 32 , 33 , 34 , 35 , 36 , 37, 38 , 39 . In presence of charge neu- 
trality between the leads and the wire and in presence of 
perfect mirror symmetry between the incoming and out- 
going leads, the conductance of odd numbered wires is 
quantized to Go, while it is less than Go for even numbered 
wires. As the charge neutrality is broken, the oscillation in 
conductance still exists, but with the distinction that for 
a given choice of charge transfer (Ac) and lead-wire hop- 
ping (i c ), the conductance values of the odd numbered 
and even numbered wires are no longer fixed quantities 
as the size of the wires is changed. For such systems, we 
further found a change of phase of conductance oscilla- 
tion from odd-even to even-odd with increasing number 
of atoms in the wire. We found that while the amplitude 
of oscillation depends on lead-wire coupling parameter t Cl 
it is the amount of charge transfer between lead and wire, 
which affects the period of oscillation. Lifting of mirror 
symmetry between two lead-wire junctions in no charge 
transfer condition is found to reduce the conductance of 
odd numbered wires below the quantized value of Go . The 
proposed technique can be easily generalized for applica- 
tion to realistic cases. To apply this approach to real sys- 
tems, one needs to generate TB parameters of lead and 
wire via self consistent calculations while multi-orbital ef- 
fect can be taken into account via multi-channel gener- 
alization of Landauer-Biittiker formula (see equation 1 17[) . 
The proposed real space technique of calculation of con- 
ductance coupled with localized Wannier basis generated 
out of self-consistent DFT calculation 41,42; can lead to a 
viable technique for study of quantum transmittance and 
conductance of nanoscale systems of various geometries in 
general. Work is in progress along this direction. 
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